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Abstract 

k> i The decay B'' — > J/ip Kg is reconstructed with J/ip — > e+e" or fj,'^fi~ and 

H ' Kg -^ TT+TT". From the fuU ALEPH dataset at LEPl of about 4 milhon hadronic 

Z decays, 23 candidates are selected with an estimated purity of 71%. They are 
used to measure the CP asymmetry of this decay, given by sin 2/3 in the Standard 
Model, with the result sin 2/3 = 0.84 ILi'q^ ± 0.16. This is combined with existing 
measurements from other experiments, and increases the confidence level that CP 
violation has been observed in this channel to 
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1 Introduction 

In the Standard Model CP violation arises from a complex phase of the quark mixing 
matrix [1[|, and this can accommodate the observed CP violation in the K sector 0. 
Precise predictions can be made of relations between asymmetries expected in B decays, 
and their detailed study will provide an important test of the model. The first step, 
however, is to establish the existence of CP violation in B decays. This has been attempted 
using inclusive methods, where the expected asymmetry is small, O{10~^), beyond the 
sensitivity of current experiments [Q]. The alternative is to use exclusive decays, where 
the asymmetries are predicted to be large, 0{1), but the branching ratios are small. 

The decay B° -^ J/ip Kg is known as the gold-plated mode for such studies, due to 
its clean experimental signature and low theoretical uncertainty. The final state is a CP 
eigenstate, to which both B° and B° can decay. The interference between their direct and 
indirect decays via B'^-B'^ mixing leads to a time- dependent CP asymmetry given by 

,,, r(BO^JMKO)-r(BO^ JMKO) . ^^ . , 

A(t) = ^—- ', Z 7=?^ , ni = -sin 2/3 sinAmat . 1 

Here r(B'' —>■ J/ipK^) represents the rate of particles that were produced as B" decaying 
to J/ip Kg at proper time t, Arrid is the oscillation frequency of the B", and (3 is an angle of 
the "unitarity triangle" of the quark mixing matrix, given by the following combination of 
matrix elements: eiTg{—VcdV*^/VtdV^l,) |^. Information can be obtained indirectly about 
the value of sin 2/?, within the context of the Standard Model, from the combination 
of other measurements that constrain the matrix elements, such as Arrid, charmless B 
decays and CP violation in the K sector. Many such fits have been made, typically 
preferring large positive values for sin 2/3 in the range 0.4-0.8 [Q; a recent example gave 
sin 2/? = 0.75 ±0.09 0. 

The first published attempt at a direct measurement was made by OPAL [[7|. They 
selected 24 candidates with an estimated purity of 60% and reported a value outside the 
physical region, sin 2/3 = 3.2 l2!o ^ O-^? which was nevertheless interpreted as favouring 
large positive values. CDF published an analysis based on 395 candidates with a purity 
of about 40%, although half of the sample has poor proper-time determination ||^. They 
measured sin 2/3 = 0.79 Il^'W (statistical and systematic errors combined). 

The key to making such a measurement at LEP is to keep the efficiency as high as 
possible. Using the latest branching ratio i3(B° -^ 3/ip K°) = (8.9 ± 1.2) x 10"^ [|], in the 
complete dataset of ALEPH, about 30 signal events are expected before the reconstruction 
efficiency is applied. The production state of the B*^ must also be determined (or 
"tagged"). The precision on sin2/? scales as 1/(1 — 2w), where w is the mistag rate, 
the fraction of incorrectly tagged events. Using a neural-network technique to combine 
the information from many observables, the lowest possible mistag rate is aimed for, whilst 
providing a tag for every event. 

In this paper, after a brief description of the ALEPH detector, the event selection is 
discussed. Details are given of the proper-time measurement, and the production-state 



tagging. The unbinned likelihood fit for the asymmetry is then presented, followed by a 
discussion of checks and systematic uncertainties. 

2 Detector 

A detailed description of the ALEPH detector can be found in P] and its performance 



m 



10 1 . Charged particles are tracked in a two-layer silicon vertex detector with double- 
sided readout (r-0 and z), surrounded by a cylindrical drift chamber and a large time 
projection chamber (TPC), together measuring up to 33 space points along the trajectory. 
These detectors are immersed in a 1.5 T axial magnetic field, providing a transverse 
momentum resolution of Ap/p = (6 x 10~^) p at high momentum (for p in GeV/c) and 
a three-dimensional impact parameter resolution of 25 /xm. The TPC also allows particle 
identification to be performed through the measurement of specific ionization (dE/dx). A 
finely segmented electromagnetic calorimeter of lead/wire-chamber sandwich construction 
surrounds the TPC. Estimators R^, Ri, and Ri are formed for electron identification, for 
the transverse and longitudinal shower shape in the calorimeter and for the dE/dx in 
the TPC, respectively; they are calculated as the difference between the measured and 
expected value for electrons, divided by the expected uncertainty. The iron return yoke 
of the magnet is instrumented with streamer tubes to form a hadron calorimeter and is 
surrounded by two additional double layers of streamer tubes to aid muon identification. 

The LEPl data were recently reprocessed using improved reconstruction algorithms. 
A new pattern recognition algorithm for the vertex detector allows groups of nearby tracks 
to be analysed together, searching for hit assignments that minimize the overall x^ of the 
event. Information from the wires and pads of the TPC are also combined to improve the 



spatial and dE/dx resolution |TT 



Monte Carlo simulated events are used to study both the signal and the background. 
The simulation is based on JETSET |T2| and is described in detail in [Q. To tune the 



selection cuts for background suppression, 6.5 million hadronic Monte Carlo events are 
used, corresponding to about 1.5 times the data statistics. In addition, a large sample 
of signal Monte Carlo events is used for the determination of the expected signal mass 
distribution, reconstruction efficiency, and for training of the neural network for tagging. 

3 Event selection 

Data taken by ALEPH in the years 1991-95 at the Z resonance are used, corresponding 
to 4.2 million hadronic Z decays. Hadronic events are selected in the data as described 
in |n. The production vertex position is reconstructed on an event-by-event basis using 



the constraint of the average beam-spot position ||T^ . 



First, the J/V' — *■ ^'''^~ reconstruction is performed. The daughter tracks are required 
to have momentum greater than 2 GeV/c, distance of closest approach to the primary 
vertex of less than 2 cm transverse to and 10 cm along the beam axis, four or more 
hits in the TPC, and polar angle satisfying | cos 6*1 < 0.95. All oppositely-charged 
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Figure 1: Reconstructed J/ipK^ mass (a) for hadronic Monte Carlo, (b) for the data; the 
superimposed fits are described in the text, and the dashed hues indicate the signal region. In 
(a) the events that truly originated from a B'^ ^ J/ipKa^ decay are shaded, and the insert shows 
the mass distribution for a large sample of signal Monte Carlo events. 



pairs of such tracks, with opening angle satisfying cos6u > 0.85, are investigated for 
lepton identification. They must both be identified as muons or electrons using loose 
identification criteria. For muons, a pattern of hits in the HCAL consistent with a muon 
is required |T^; for electrons, cuts are made on the estimators: |i?T| < 4, i?L > —3 
and Ri > —4. The invariant mass of the lepton pair is required to be in the range 2.6- 
3.3 GeV/c^, with the large window around the J/ip mass (particularly on the low side) 
maintaining high efficiency in the presence of radiative decays or bremsstrahlung. 

Next the Kg -^ tt'^tt" reconstruction is performed, as described in |T^. The distance 
of closest approach of the two daughters must be less than 5 mm, and the dE/dx for each 
daughter is required to be within three standard deviations of the expected value for a 
pion. The angle between the reconstructed directions of the Kg and ^/ip is required to 
satisfy cos6'^k > 0.85, and the resultant of their momenta must satisfy pb > 22 GeV/c. 
The reconstructed Kg invariant mass is required to be within 15 MeV/c^ of its nominal 
value. 

A fit is made for the B° decay vertex, using the two lepton tracks and the Kg. The 
vertex is required to be successfully reconstructed, and the decay length is determined from 
the distance between the production and decay vertices, projected along the momentum 
vector of the B° candidate. The decay length is required to be greater than —1mm, 
and its calculated error less than 1 mm. Finally, the J/V' Kg invariant mass is calculated, 
assigning the nominal masses to the two particles. The resulting distribution is shown 
as an insert in Fig. |1| (a) for signal Monte Carlo. The peak is fitted with the sum of two 
Gaussian functions, the first accounting for 72% of the events with a width of 20 MeV/c^ 
and the second with width 90MeV/c^, offset to lower mass by 50MeV/c^. The signal 
region is defined as 5.05-5.45 GeV/c^, and the overall reconstruction efficiency is 28%. 

When the event selection is applied to the hadronic Monte Carlo sample, the resulting 



J/V' Kg mass distribution is shown in Fig ^(a). There are 20 events in the signal region, of 
which 11 are from background. The superimposed fit is the sum of an exponential shape to 
describe the background and the signal shape discussed above. The four fitted parameters 
are the background slope and normalization, and the signal mass and normalization. 

The mass distribution that results when the event selection is applied to the data 
is shown in Fig. |l](b). A clear signal is seen for the B", with 23 events in the signal 
region. A fit similar to that of the hadronic Monte Carlo is made, giving a fitted B° 
mass of 5.26 ± 0.01 GeV/c^, slightly lower than, but consistent with, the world-average 
value of me = 5.279 GeV/c^. The fit is used to assign an event-by-event background 
probability for each event in the signal region, which is used in the maximum likelihood 
fit for the CP asymmetry. Their sum corresponds to 6.6 background events, giving an 
average background fraction /bkg = 0.29 ± 0.06, where the uncertainty is estimated by 
varying the parametrizations within their statistical errors, and using alternative shapes 
for the background. The fitted shape and normalization of the background are consistent 
with those seen in the Monte Carlo. 

After subtraction of the background, about 16 signal events remain, to be compared 
with the predicted number of 9 ± 2, calculated from the expected production rate and 
decay branching ratios. Figure ^ shows a particularly clean signal candidate, in which 
there are no other charged tracks in the signal hemisphere. 

4 Proper-time determination 

The decay length of the B° is determined from vertex reconstruction as described in 
the previous section. The decay-length resolution determined using the Monte Carlo 
simulation is reasonably described by a Gaussian distribution of width 200 /xm, although 
there are small tails from events with poorly measured vertices. The uncertainty on the 
decay length is estimated by propagating the production and decay vertex errors, which 
are calculated in turn from the errors of the tracks used in their determination. The pull 
distribution, given by the difference of true and reconstructed decay lengths divided by 
the calculated decay-length uncertainty, is close to being normally distributed. It has a 
fitted Gaussian width of 1.1, indicating that the event-by-event estimate of decay-length 
uncertainty is reasonably accurate. 

As the signal events are fully reconstructed, the momentum resolution is excellent. 
There is a Gaussian core of 1% relative error, with an overall RMS of 2.5%. This small 
uncertainty on the momentum measurement gives a significant contribution to the overall 
proper-time uncertainty only for events with long decay lengths, greater than about 1 cm. 

The proper time is given by t = d-Qm-Q/p-Q. Its uncertainty is calculated as follows: 




^t = t\r-^\ +^ . (2) 



The first term in parentheses is taken from the event-by-event measurement of the decay 
length d-Q and its error cr^, scaled by a factor 1.1 ± 0.2, where the uncertainty is taken 
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Figure 2: Event display of a B° ^ 3/ipK^ candidate in the data, with J/^ -^ e+e": (a) 
fish-eye view of the projection transverse to the beam axis, (b) zoom into the vertex region for 
the same projection, with reconstructed vertices marked with elhpses, and the reconstructed 
neutrals marked with dashed lines. 



into account for systematic studies; the second term is taken as (2 ± 2)%. Despite being 
conservative, these estimates of the uncertainty lead to a negligible effect on the measured 
asymmetry, as the characteristic scale of the proper-time development of the asymmetry 
is l/Ar?T,d ~ 2ps, much longer than the typical proper-time resolution of 0.1 ps. 

Of the 11 background events in the signal region for the hadronic Monte Carlo, 
nine involve tracks from b-hadron decays. The apparent proper time and its error are 
determined for these background events, following the procedure discussed above, and a 
fit made for their effective lifetime Tbkg. The probability density function is taken as the 
convolution of an exponential lifetime distribution with a Gaussian resolution function, 
and the fit gives Tbkg = 1.4 toJps. 



5 Production-state tagging 

The extraction of the CP asymmetry from the signal events requires the determination of 
whether they originated from a B° or B° at production. This is achieved by studying the 
properties of the rest of the event, excluding the lepton and pion pairs that come from the 
signal B*^ decay. For this purpose, two hemispheres are defined with respect to the thrust 
axis of the event, determined using both charged and neutral energy-flow objects |I0| . 
These are used to separate information from the same and opposite sides of the event 
with respect to the B*^ candidate. Properties of the opposite side are used to determine 
the particle/antiparticle nature of the other b hadron that was produced in conjunction 
with the signal B°, and thus indirectly determine its production state. The same side 
carries information from the fragmentation process that produced the signal B°, which 
can also be used in the production-state determination. 

The calculation of the tag for the opposite side starts with the search for a secondary 
vertex due to b-hadron decay. This is achieved using a topological vertexing algorithm that 
combines information from all charged tracks in the hemisphere. Jets are reconstructed 



from the charged tracks and neutral energy-flow objects using the JADE algorithm ||17| , 
with a jet-resolution parameter of 0.02. The highest energy jet and the jet which forms 
the highest invariant mass with it are selected, and the secondary vertex is constrained 
to lie along the direction of the selected jet in its hemisphere, within errors. Each track 
is then assigned a relative probability V^ that it originates from the secondary vertex. 

Next, the b-hadron flight direction is estimated. Jets are reconstructed with a jet- 
resolution parameter lowered to 0.0044, as this gives an improved estimate of the b-hadron 



direction [1^. If more than one jet is found in the hemisphere considered, the b-jet 
candidate is chosen on the basis of the kinematic properties of its tracks and the presence 
of lepton candidates. The leading track of the b-jet candidate is then used as a seed for 



a second cone-based jet algorithm |]T8|, which is taken as a first estimate of the b-hadron 
flight direction. In two-jet events, the thrust axis is chosen instead. The uncertainties on 
the reconstructed angles are parametrized from the simulation, as a function of the jet 
momentum. A second estimate of the direction is taken as the vector joining the primary 
and secondary vertices, and its error is parametrized as a function of the measured decay 
length. The two estimates are averaged using their parametrized errors, and the result 
taken as the b-hadron flight direction. 

To construct the charge estimators, b-decay and fragmentation tracks must be 
distinguished. This is achieved by combining three variables for each track in the 
hemisphere: the rapidity and longitudinal momentum relative to the b-hadron flight 
direction, and the vertex assignment probability V^. They are combined in a neural 
network, along with a hemisphere-based indicator of the quality of the inclusive vertex 
finding, which acts as a "control variable". This does not directly discriminate between 
the b-decay and fragmentation tracks, but improves the overall performance of the neural 
network as the discriminating power of the other inputs vary as a function of the control 
variable. The single track separation output is shown in Fig. |^(a), using neural networks 
trained separately for classes of tracks passing different quality criteria. The neural- 



network output value Xsec is converted to a probability Psec that the track comes from the 
secondary vertex. The b-hadron momentum is then determined from the sum of charged 
energy assigned to the secondary vertex (based on Psec weights), the longitudinal fraction 
of neutral energy in the jet, and the missing energy in the hemisphere. 

Nine charge- sensitive inputs are used for the opposite-side production-state tag: 

1. Jet charge (k, = 0.5): 

^'- E.iPrd,r ' ^'^ 

where i runs over the charged tracks in the hemisphere, g, and Pi are the charge 
and momentum of the track, dj is the b-hadron flight direction. 

2. Jet charge (k, = l.Oj.' defined as in Eq. 3. 

3. Charge of the track with highest longitudinal momentum in the hemisphere, 
calculated relative to the b-hadron flight direction. 

4. Total charge: Qtot = J2i li where i runs over the charged tracks in the b jet. 

5. Secondary vertex charge: Qvtx = Z^i "^sec ^i where i runs over all charged tracks in 
the hemisphere. 

6. Weighted primary vertex charge: 

^"'" j:.ii-vuiPrdjr ' ^ ^ 

where i runs over all charged tracks in the hemisphere and k, = 1.0. 

7. Weighted secondary vertex charge: as for Qpri but replacing (1 — Pgec) with V^^^, and 
taking k, = 0.3. 

8. Decay kaon charge: kaons are identified using dE/dx, based upon the ratio of 
probabilities that the measured ionization is due to a kaon, relative to either a 
kaon or a pion. This variable is combined with the track momentum, longitudinal 
momentum and Psec value to select kaons from the b — ;> c — *> s cascade using a 
further neural network. The reconstructed b-hadron momentum in the hemisphere 
is also used as a control variable. The charge of the track with the highest output 
from the neural network is used to sign the output value. 

9. Lepton charge: lepton identification is performed with tighter requirements than 
those described previously [^ . If more than one lepton candidate is selected, that 



with the highest p^ with respect to the jet axis is used. The lepton transverse 
momentum is signed by the track charge. 
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Figure 3: Neural-network output distributions for signal Monte Carlo events: (a) separation 
of b-decay and fragmentation tracks, Xseo with the contribution from b-decay tracks shaded; 
(b) opposite-side production state, Xopp, with the contribution from hemispheres containing b 
hadrons shaded; (c) overall event production state, 2;tag, with the contribution from B" decays 
shaded, (d) Calibration of the production-state tag q versus the neural-network output, with 
superimposed linear fit; the effect of degrading the mistag rate for systematic studies is indicated 
by the dashed line. 



Finally, six additional control variables are used: 
multiplicity, the spread in the measured values of "Pse. 



cos 6'thrust, the charged track 



o-s 



L-ii ' sec \^ 'SI 



the 



reconstructed b-hadron momentum, the reconstructed decay length, and the lepton 
momentum (if a lepton has been selected). These are combined along with the charge 
estimators using a neural network, which takes into account correlations between the 
variables. It provides a single output, Xopp, which is shown in Fig. ||(b). 

Information on the production state from the same hemisphere as the B*^ candidate 
is limited to the charged tracks from fragmentation. Excluding the tracks from the B°, 
seven charge estimators are constructed using tracks coming from the primary vertex in 
the signal hemisphere. In a similar way to the opposite-side analysis, these include two jet 
charges with k values of 0.5 and 1.0, the charge of the track with the highest longitudinal 
momentum with respect to the jet, one momentum-weighted and two longitudinal- 
momentum- weighted primary vertex charges with k values of 1.0, 0.3 and 1.0, respectively. 



and finally the sum of track charges inside the jet. The output of the opposite-side neural 
network, Xopp, is then combined with these same-side charge estimators and the following 
control variables: (Jsec, cos6'thrust, and the same-side charged track multiplicity, using 
another neural network. The output of this event tag, Xtag, is shown in Fig. 0(c). 

If events with Xtag > 0.5 are taken to have a B° in the production state and those with 
^^tag < 0.5 to have a B°, the fraction of incorrect tags (the average mistag rate) is 27.8%. 
The event-by-event value of the neural-network output a;tag is used in the measurement of 
the asymmetry. Due to its peaked distribution, seen in Fig. ^(c), the effective mistag rate 
is lower than the average value quoted above. Using an overlap integral technique ||T9| 



the effective mistag rate is found to be 24.1%, measured with an independent sample of 
simulated events, with an expected difference of (0.2 ±0.2)% between the effective mistag 
values for B° and B° hemispheres. The opposite-side tag alone gives average and effective 
mistag rates of 31.4% and 27.6%, respectively. 

Finally, the production-state tag q is calculated from the neural- network output Xtag, 
correcting for the purity in each bin of Fig. ^(c) using 

^ F(xtag)+G(Xtag) ' ^^ 

where F(xtag) is the distribution of a;tag for events produced as B°, shaded in the figure, 
and G{xta.g) is the distribution for B"; q = +1 for events produced as B*^ and —1 for B" in 
the case of perfect tagging. The relationship between Xtag and q, determined with signal 
Monte Carlo events, is consistent with linearity, as shown in Fig. ^(d). It is parametrized 
as a straight line q = atag + 2btag{xtag — 0.5), with fitted coefficients atag = —0.01 ± 0.01 
andfetag = 0.95 ±0.01. 

6 Asymmetry measurement 

The production-state tag q is plotted against the proper time t in Fig. ||(a), for the 23 
events in the signal region of the data. Considering the events which have a clear tag 
result, |g| > 0.5, there are 4 "B°-like" events (positive q) and 9 "B°-like". Furthermore, 
the excess of events with negative q is more noticeable with increasing proper time, as 
would be expected for a negative CP asymmetry. 

To extract a measurement of the asymmetry, an unbinned likelihood is calculated. 
The probability density function expected for the signal distribution is given by 

g-t/f"d 

'Psigit, q) = — (1 - q sin 2(3 sin Arudt) . (6) 

2rd 

The B° lifetime Td = 1.548±0.032ps and Amj = 0.472±0.017ps-^ are fixed to the central 
values of their world averages [Q (the uncertainties are taken into account in systematic 
error studies). 

A convolution is made of this signal distribution and a Gaussian resolution function 
TZ{at) with width given by the event-by-event proper-time resolution calculated in 
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Figure 4: (a) Production-state tag q versus proper time t for the 23 events in tlie signal region 
of the data. The proper-time uncertainty is indicated for each point, and histograms of the two 
variables are shaded. The solid points indicate events with a background probability less than 
30%. (b) Expected proper-time distribution of B*^ and B'' events, for fixed values of the CP 
asymmetry and proper-time resolution: sin 2/? = 0.7, at = 0.2 ps. 

Section 4. The result is illustrated in Fig. §(b), where the contributions from B'' and 
B° are indicated as a function of proper time for fixed values of the asymmetry and 
resolution. 

The probability density function for background events Vhkg{t, q) is taken to have the 
same form as that of the signal but replacing t^ -^ Tbkg and sin 2/5 -^ abkg, where the 
effective lifetime of the background Xbkg was determined in Section 4. The background 
asymmetry Obkg is taken to be zero, but is varied to study possible systematic effects. The 
likelihood of an event is then calculated as 



jC-i = (1 - /bkg) Vsigit, q) ® n{at) + /bkg ^bkg(i, q) ® TlicTt) , 



(7) 



where /bkg is the event-by-event background fraction discussed in Section 3. A generalized 
likelihood function is used, adding a Poisson term to the combined likelihoods of the signal 
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Figure 5: (a) Log-likelihood versus sin 2/3 from the fit, with respect to the minimum, 
(b) Reconstructed asymmetry versus proper time, after correction for the average dilution in 
each time bin. The result of the fit for the asymmetry is shown superimposed; note that this is 
not fitted directly to the binned data shown. 



candidates: 



u e~ 



n\ 



HA 



0.84l?:^2^ 



where n = 23 is the number of observed candidates, and v is the number expected. 

A fit is made with two free parameters, sin 2/5 and z/, with the results sin 2/5 
V = 23.3 lll'l. A scan of the log-likelihood versus sin 2/5 is shown in Fig. |^(a). The result 
of the measurement is displayed as a function of proper time in Fig. ^(b). For the purpose 
of display the data are divided into three proper-time bins, and the average asymmetry is 
corrected for the average dilution (from mistagging and background) in each bin. Negative 
asymmetry, and thus positive sin 2/5, is favoured. 



7 Checks and systematic uncertainties 

The statistical error from the fit has been checked using a fast Monte Carlo simulation. 
Event samples the same size as the data were generated according to Eq. 6 and 7, using the 
resolution function and tagging distributions discussed above. The value of sin 2/5 was set 
to that measured for the data, and all other parameters were fixed to their nominal values 
in the fit. The number of background events within the sample was varied according to 
Poisson statistics. The CP asymmetry of each sample was measured as if it were data, 
and the central value and errors recorded. This was repeated for many samples. The 
values reconstructed for the positive and negative errors of the data are consistent with 
the distributions of values seen in the simulated samples. Furthermore, the measured 
errors reasonably estimate the spread of the measured central value about its true input 
value: the pull is normally distributed. 

The contributions to the systematic error are listed in Table 1. The lifetime and 
oscillation frequency of the signal were varied within their world-average uncertainties. 
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Table 1: Contributions to the systematic uncertainty on sin 2/3. 



Source 


Variation 


'^sys 


B° lifetime 


rd = 1.548 ± 0.032 ps 


<0.01 


B° oscillation 


Amd = 0.472 ±0.017ps-i 


0.04 


Decay-length resolution 


aa X (1.1 ±0.2) 


<0.01 


Momentum resolution 


Ap/p = 0.02 ± 0.02 


<0.01 


Background level 


/bkg = 0.29 ± 0.06 


0.09 


Background lifetime 


rbkg = 1.4l|]:ips 


0.01 


Background asymmetry 


abkg = 0.0 ± 0.2 


0.04 


Tag calibration offset 


fltag = -0.01 ± 0.02 


0.10 


Tag calibration slope 


6tag = 0.95 ± 0.05 


0.05 


Total 




0.16 



and the effect on the measured asymmetry taken as a systematic error. For the lifetime, 
a check was made by fitting for the lifetime of selected signal Monte Carlo events, giving 
r = 1.56 ± 0.02 ps, in agreement with the input value of 1.56 ps, indicating no evidence 
for a systematic bias. Similarly their asymmetry was measured to be —0.02 ± 0.03, in 
agreement with the input value of zero. The decay-length and momentum resolution 
parametrizations were varied as discussed in Section 4. 

For the background, the uncertainty on the level was discussed in Section 3. The 
effective lifetime of the background was varied within the large uncertainty measured 
using the hadronic Monte Carlo, given in Section 4. To increase the statistics for 
the study of the CP asymmetry of the background, events in the side-band region 
4.0 < m(J/'?/'Kg) < 5.0GeV/c^ were investigated. There are 140 such events in the 
data, with asymmetry 0.09 ± 0.23, indicating no significant effect. About one event is 
expected in the signal region from the decay B° — > J ftp K*^ with K*" —>■ Kg tt^, predicted 
to be mainly CP even. This would correspond to a background asymmetry of less than 
0.1. For the systematic error evaluation the asymmetry of the background, flbkg, was 
varied by ±0.2. 

Uncertainties arising from the production-state tag can be considered as being due to 
knowledge of the overall mistag rate for signal events in data and the possible difference in 
the individual mistag rates for B° and B'' events. Samples of events were isolated in data 
with a similar topology to that of the signal, where there is a b hadron of known charge 
with decay products that can be cleanly separated from fragmentation tracks. Three 
samples of B^ candidates were selected for this purpose. The same neural network was 
used for tagging their production state as for B° — > J/^Kg, except that the sign of the 
same-side charge estimators is reversed. This is required as a b quark combines with a d 
(u) quark to produce a B° (B+), and the accompanying d (u) quarks have opposite sign. 

1. B+ — » J/ipK'^ decays: The signal selection is modified slightly, to select B"*" — > 
J/ipK'^ decays (and charge conjugate). The reconstructed J/t/^K"*" mass in the data 
is shown in Fig. ^(a); a clear signal is seen, with 52 events in the signal region 
compared to an expected background of about 17. In this channel the efficiency 
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Figure 6: (a) Reconstructed J/^ K+ mass, for the data, (b) Production-state tag neural-network 
output Xtag, when applied to events from the signal region. 

is similar to that of J/?/' Kg, but the product branching ratio is about three times 
higher. Now the charge of the kaon indicates whether the signal was produced by 
B"*" or B~, and the resulting distributions of the neural- network output Xtag are 
compared in Fig. §|(b). From a study of side-band events, the background is found 
to have a similar tagging distribution to the signal. The average mistag rate is 26% 
in signal Monte Carlo and is measured to be (27 ±6)% in the data, giving confidence 
that the Monte Carlo simulation reproduces faithfully the features of the data used 
for tagging. The asymmetry and lifetime have also been measured for these events, 
with the results 0.09 ± 0.41 and 1.57 io'2oPS, respectively. The latter agrees with 
the world-average value of 1.653 ± 0.028 ps p|]. 

2. B+ -^ VPt^X decays: A selection is made of D° candidates, decaying to a kaon and 
a pion of opposite charge. Only those D° candidates are used which are found to lie 
in the same hemisphere as an identified lepton (e or /i) with transverse momentum 
greater than 0.5 GeV/c and charge opposite in sign to that of the pion from the D" 
decay. The estimated b purity of the sample is 87%, dominated by B='= decays. 
The lepton charge indicates whether the signal was produced by a B+ or B~, 
and the resulting tag distributions are shown in Fig. 0(a). The effective mistag 
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Figure 7: Production-state tag neural-network output for (a) D''^"'" and (b) vertex-charge 
selected samples. In each figure the upper plot is for events selected as B"*" candidates, the 
lower for B~. 



rate determined from data is (24.5 ± 0. 
(24.1 ±0.5)%. 



3. 



with a Monte Carlo expectation of 



Inclusive B^ decays using vertex charge: A selection of B^ decays is made by 
requiring a hemisphere with o"sec < 0.32 opposite to a b-tagged hemisphere pO| , 
that gives an estimated b purity of 95%. A cut is then made on the absolute value 
of the vertex charge |Qvtx| > 0.6 in the selected hemisphere, which isolates a sample 
containing an expected 73% of B^ decays. The sign of the vertex charge indicates 
whether the signal is produced by a B^ or a B^ and the resulting distributions of 
the tag output are shown in Fig. 0(b). Although significant differences are seen 
in this sample between the tag distributions for B"*" and B~ candidates, due to 
material effects, they are well reproduced by the Monte Carlo. The effective mistag 
rate determined from data is (21.9 ± 0.2)%, with a Monte Carlo expectation of 
(21.0 ±0.1 



70. 



The largest data-Monte Carlo disagreement for the mistag rates is (0.9 ± 0.2)% from 
the vertex-charge selected sample. In addition, separating the samples into B"*" and B~ 
decays as shown in Fig. and determining the mistag rates separately for each results in a 
maximum observed discrepancy of (2.4 ± 1.3)% from the D^^"*" selected events. These are 
taken as systematic uncertainties for the production-state mistag rate of the B° -^ J/ip Kg 
signal. Their effect is propagated to the fit for the CP asymmetry by modifying the 
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Figure 8: Likelihood distributions versus sin 2/? for the three existing measurements (dashed 
lines) and their combination (shaded). 



calibration of the tag, the parametrization of q versus Xtag- A degraded mistag rate leads 
to a reduced slope of the calibration, as illustrated in Fig. ^(d). Similarly, a B°-B° 
tagging difference is seen as an offset to the calibration. The variations applied to the 
parametrization are listed in Table 1. 

Adding the systematic error contributions in quadrature, the final measurement is 



sin2/5 = 0.84+5:^2 ^0.16 , 
where the first error is statistical and the second systematic. 



(9) 



8 Conclusion 



An analysis of B° -^ J/ip Kg decays has been performed, with J/ip — > e"'"e~ or /x+/x^ and 
Kg -^ TT+vr^. A reconstruction efficiency of 28% is achieved. From the full dataset of 
ALEPH at LEPl of 4.2 million hadronic Z decays, 23 candidates are selected with an 
estimated purity of 71%. They are used to measure the CP asymmetry of this decay, 
with the result sin 2/3 = 0.84 t^'xtt, where the uncertainty is dominated by the statistics. 
This result is compared with the other two published measurements in Fig. §. The 
likelihood of each measurement is approximated as a Gaussian distribution on either 
side of the central value, with width equal to the sum in quadrature of statistical and 
systematic errors. The three log-likelihoods can be summed to give a combined result, 
i.e. neglecting any correlation between the systematic errors of the different experiments 
(expected to be small). The resulting likelihood is shown by the shaded distribution in the 
figure and corresponds to sin 2/5 = 0.88 ^^39. The integral of the likelihood for sin 2/3 > 
is 77% for the ALEPH result alone, or 67% if the total integral is limited to the physical 
region |sin2/3| < 1. The corresponding fractions are 98.6% and 98.0% for the combination 
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of the three analyses. Thus the confidence level that CP violation has been seen in this 
channel is increased to 98%, compared to 93% for CDF alone P|. 



Preliminary results from the BABAR and BELLE experiments |21| are consistent with 



the combined result given above, with similar precision and slightly lower central values. 

Acknowledgements 

It is a pleasure to thank the CERN accelerator divisions for the successful operation of 
LEP. We are also grateful to the engineers and technicians in our institutes for their 
contribution to the performance of ALEPH. Those of us from non-member states thank 
CERN for its hospitality. 

References 

[1] M. Kobayashi and K. Maskawa, Prog. Theor. Phys. 49 (1973) 652. 

[2] Particle Data Group, D. Groom et a/., Eur. Phys. J. C 15 (2000) L 

[3] OPAL Collaboration, K. Ackerstaff et ai, "A study of B meson oscillations using 
inclusive lepton events in hadronic Z decays", Z. Phys. C 76 (1997) 401; OPAL 
Collaboration, G. Abbiendi et ai, "Measurement of the B"*" and B° lifetimes and 
search for CP(T) violation using reconstructed secondary vertices", Eur. Phys. J. 
C 12 (2000) 609; DELPHI Collaboration, M. Feindt et al, "Study of inclusive 
CP(T) asymmetries in B decays", DELPHI 97-98 CONE 80, contributed paper 
449 to HEP97 Jerusalem (August 1997); ALEPH Collaboration, R. Barate et a/., 
"Investigation of inclusive CP asymmetries in B° decays" , submitted to Eur. Phys. 
J. (July 2000). 

[4] See, for example: Y. Nir and H. Quinn, in "B Decays", ed. S. Stone, World Scientific 
(1994) 362; C. Dib, I. Dunietz, F. Oilman and Y. Nir, Phys. Rev. D 41 (1990) 1522. 

[5] See, for example: A. Ali and B. Kayser, "Quark mixing and CP violation" , in "The 
Particle Century", ed. Gordon Eraser, Inst, of Physics Publ. (1998) |hep-ph/9806230| ; 



F. Caravaglios, F. Parodi, P. Roudeau and A. Stocchi, "Determination of the CKM 
unitarity triangle parameters by end 1999", to appear in Proc. of 3rd Int. Conf. 
on B Physics and CP Violation (BCONF99), Taipei, Taiwan (December 1999) [hep 
ph/0002171| . 



[6] S. Mele, Phys. Rev. D 59 (1999) 113011. 

[7] OPAL Collaboration, K. Ackerstaff et ai, "Investigation of CP violation in B*^ -^ 
J/V'Kg decays at LEP", Eur. Phys. J. C 5 (1998) 379. 

[8] CDF Collaboration, T. Affolder et al, "A Measurement of sin 2/5 from B -^ J/^Kg 
with the CDF Detector", Phys. Rev. D 61 (2000) 072005. 

16 



[9] ALEPH Collaboration, D. Decamp et ai, "ALEPH: a detector for electron-positron 
annihilations at LEP", Nucl. Instr. Methods A 294 (1990) 121; B. Mours et al, 
Nucl. Instr. Methods A 379 (1996) 101. 

[10] ALEPH Collaboration, D. Buskulic et ai, "Performance of the ALEPH detector at 
LEP", Nucl. Instr. Methods A 360 (1995) 481. 

[11] ALEPH Collaboration, R. Barate et ai, "Measurement of the B'' and B+ meson 
lifetimes", CERN-EP/2000-106, submitted to Phys. Lett. (August 2000). 

[12] T. Sjostrand and M. Bengtsson, Comp. Phys. Comm. 43 (1987) 367. 

[13] ALEPH Collaboration, D. Buskulic et ai, "Heavy flavour production and decay with 
prompt leptons in the ALEPH detector", Z. Phys. C 62 (1994) 179. 

[14] ALEPH Collaboration, R. Barate et ai, "Measurement of the Z resonance parameters 
at LEP", Eur. Phys. J. C 14 (2000) 1. 

[15] ALEPH Collaboration, D. Buskulic et ai, "Heavy quark tagging with leptons in the 
ALEPH detector", Nucl. Instr. Methods A 346 (1994) 461. 

[16] ALEPH Collaboration, D. Buskulic et ai, "Production of K° and A in hadronic Z 
decays", Z. Phys. C 64 (1994) 361. 

[17] JADE Collaboration, W. Bartel et ai, "Experimental studies on multijet production 
in e+e" annihilation at PETRA energies", Z. Phys. C 33 (1986) 23. 

[18] OPAL Collaboration, R. Akers et ai, "QCD studies using a cone-based jet finding 
algorithm for e+e" collisions at LEP", Z. Phys. C 63 (1994) 197. 

[19] ALEPH Collaboration, D. Buskulic et ai, "Study of the B^-B^ oscillation frequency 
using D-£+ combinations in Z decays", Phys. Lett. B 377 (1996) 205; D. Jaffe, F. Le 
Diberder and M.-H. Schune, LAL 94-67 FSU-SCRI 94-101. 

[20] ALEPH Collaboration, R. Barate et ai, "Search for the Standard Model Higgs boson 
at the LEP2 collider near v^ = 183 GeV", Phys. Lett. B 447 (1999) 336. 

[21] BABAR Collaboration, B. Aubert et ai, "A study of time dependent CP violating 
asymmetries in B° -^ J/V'Kg and B° -^ V(2S) K° decays", SLAC-PUB-8540, |eF 



ex/0008048| (July 2000); H. Aihara, "New results from BELLE", presented at XXX 
Int. Conf. on HEP, Osaka, Japan (July 2000). 



17 



